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ABSTRACT 


The problem of an elastic solid containing a semi-infi- 
nite plane crack subjected to concentrated shears parallel 
to the edge of the crack is considered in this paper. A 
closed form solution using four distinct harmonic functions 
(none of which can be taken arbitrarily) is found to satisfy 
the finite displacement and inverse square root stress singu- 
larity at the edge of the crack. Explicit expressions in 
terms of elementary functions are given for the distribution 
of stress and displacement in the solid. These are obtained 
by employing Fourier and Kontorovi ch-Lebedev integral trans- 
forms and certain singular solutions of Laplace equations in 
three dimensions. The variations of the intensity of the 
local stress field along the crack border are shown graphical 
ly. The present analysis offers an example which is in con- 
trast with the conclusion established in the literature that 
one of the four Papkovi ch-Neuber functions in three-dimension 
al elasticity may be arbitrarily set to zero. 
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INTRODUCTION 


The general notion of the classical theory of elasticity 
* 

[1-3] is that one of the four harmonic functions appearing 
in the Papkovi ch-Neuber solution of the displacement equations 
of equilibrium may be set to zero (or arbitrarily chosen) 
without loss in completeness of the solution. A counter ex- 
ample is provided in this paper on an unfinished problem ini- 
tiated by Uflyand [4]. The problem, basically, is concerned 
with determining the distribution of stress in a half-space 
when its plane surface is subjected to mixed conditions sepa- 
rated by an infinite rectilinear boundary. 

The solution to the problem of concentrated normal and 
shear forces applied to the surfaces of a semi -i nf i ni te plane 
crack serves as the Green's function for a number of crack 
problems of interest in fracture mechanics. A sketch of the 
crack configuration is shown in Figure (1). In terms of the 
coordinates indicated, the surfaces of the crack are defined 
by 6 = +tt , 0<r<°°, -co<z<co. Without affecting the generality 
of the problem, the points of application of the loads may be 
taken as r = a, 0 = +ir, z = 0. Because of skew-symmetry, it 
suffices to consider a half-space y>_0 with mixed conditions 
specified on the plane y = 0. It was Uflyand [4] who made 
the first attempt to solve this problem by using only three 
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harmonic functions in the space y>0 and the Kontorovich- 
Lebedev integral transform in the variable r. The fourth 
function in the classical Pap kovi ch-Neuber representation of 
the displacement field was arbitrarily neglected. This tech- 
nique has failed to yield a solution for the case of two equal 
and opposite force R directed parallel to the edge of the 
crack [4, page 382]. 

The objective of this paper is two fold. First, the 
Green's function for the three-dimensional crack problem is 
constructed. Next, it is shown that the fourth Papkovich- 
Neuber function contributes to the solution and is not zero. 
Four distinct harmonic functions are required to satisfy the 
finite displacement and the inverse square root of r condi- 
tions as r-*0. The state of displacement and stress throughout 
the solid is determined in closed form in terms of elementary 
functions. These are obtained by employing a Fourier trans- 
form in the variable z, Kontorovich-Lebedev transform in r 
and singular solutions of the Laplace equation in three dimen- 
sions [5]. The stress -i ntens i ty factors due to the loading R 
are determined explicitly and their variations along the crack 
border are shown graphically. 

BASIC EQUATIONS 

The equilibrium of an isotropic and homogeneous solid is 
governed by the Navier displacement equations which in the 
absence of body forces and in vector notations appear as 
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( 1 ) 


V 2 u + VV. u = 0 

where u is the displacement vector, V is the gradient operator 
and v designates Poisson's ratio of the material. Denoting 
the projections of the displacement in the directions of cy- 
lindrical coordinates by ( u ^ , u g , u z ), the Papkovich-Neuber 
general representation of equations (1) gives the displacement 
field 


= 4(l-v)(f-|COse + f 2 s i n 0 ) 

3 F 
3r 

(2a) 

= 4(1 -v) (f ^cose - f-jSine) 

1 3 F 
’ r 36 

(2b) 

■ 4(l-v)f 3 - |f 


(2c) 


in which the following abbreviation has been introduced 

F = fg + (rcosG)f-] + (rsine)f£ + zfg (3) 

and f (r,0,z), n = 0,1, 2, 3, are space harmonic functions. The 
corresponding stress field is readily obtained from equations 
(2) and (3) and the usual stress-displacement relations in 
linear el astostatics. In particular the stresses associated 
with the 0-plane are found: 
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a fl 3f l df 9 

2 ^J = - (1 -2v) (cose Jf- + Sine yp- 


^r’ V) (cos6 a/ ‘ s1n0 ¥9 


3 f , 3 f _ 

’■> + 2v 


82f ( 


3 2f 0 1 3 2 f -i 

- - (cose 3 pr + sine 


3 2 f 


3 z ' 


36 
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3 2 f . 


3 2 f 


+ z ( 3r ^ + 3 i 
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(4a) 


re 


2y 


li + 

3 r 


Hlzvl (cos0 __1 + sin0 __2 
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(4b) 


6z 


2y 


li . 2 ( 1 -y ) 3 f 3 

3 z r 36 


(4c) 


In equations (4), y depicts the shear modulus of the solid and 
the function G(r,e,z) is defined through the relation 


G = ( 1 -2v) (f ^cose - f^sine) - 


1 11° 

r 30 


3 f 


- COS 


1 


30 


- sine 


3f , 

t 

36 


zOl 

r 36 


(5) 


The remaining stress components a z , and r r2 are not needed 
for the mere purpose of completing the analysis and will not 
be mentioned. 

THE PLANE CRACK PROBLEM 


Consider the case of a pair of equal and opposite concen- 
trated shear forces R applied to the surfaces of a plane as 
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shown in Figure 1 (with P = Q = 0). For this problem, the 
displacements and stresses are skew-symmetric with respect to 
the variable z. Hence, u r , u , a 0 and x rQ are odd in z while 
u and t Q are even in the same variable. In addition, the 
deformed solid also exhibits skew-symmetry in the coordinate 
0, namely u r , u 2 , a Q being odd in 0 and u Q , x r0 , r Qz even in 
0. The latter condition suggests that the problem can be 
formulated for the upper half-space y>0 with appropriate con- 
ditions prescribed on the planes 0 = 0 and 0 = tt. In observ- 
ing these conditions, the continuity of the solid across the 
plane 0 = 0 implies 


u r (r,o,z) = 0 

(6a ) 

u z (r,o,z) = 0 

(6b) 

Q 

CD 

-5 

O 

W 

N 

II 

O 

(6c) 

On the upper surface of the crack 0 = tt, 
scribed by 

the loading is d e - 

x 0z (r,TT,z) = T 0 (r,z) 

(7a) 

T r0 (r,TT,z) = 0 

(7b) 

( r ,tt , z ) = 0 

(7c) 
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in which T Q (r,z) is a specified function describing the ap- 
plied shear loading. Aside from the requirements given by 
equations (6) and (7), the regularity conditions at infinity 

must also be satisfied, i.e., the displacements and stresses 

- 1 - 2 

must behave as L and L when 
L = (r 2 +z 2 ) 1/2 - - 


Furthermore, near the crack boundary (r->-0), the displacements 

must be finite and the stresses are expected to have the usual 

- 1/2 

square root singularity, r ' . 


Inserting relations (6) into equations (2a), (2c) and 
(4a) and making use of (3), the functions f (n = 0,1, 2, 3) 
can be determined from the following system of equations in 
the region 0=0 



9f o 

9 r 

3f i 
r 3r 

9 f o 

( 3-4v ) f i - 



9 f 0 

9 z 

9f l 
r 9z 

9f , 
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- Z 3Z 3 
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0 


(l-2v) 


3f oti i 3 f o 
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9r r 90 


2v 


9f . 

9Z~ 


3*f ( 
9 r 2 


° f ( 

9^ 


1 32f l 

r 90 2 “ 


9 2 f Q 9 2 fo 
z( 9~P~ + lyz 7 ^ = 0 


(8a) 

(8b) 


(8c) 


Upon recognizing the identity 
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(9) 


1 32f l 3f l 


r 36' 3r 


+ + 0 = 0 


equations (8) are satisfied by requiring 


f g ( r ,0 , z ) = 

0 

(10a) 

f 1 ( r ,0 ,z ) = 

0 

(10b) 

f 3 ( r ,o , z ) = 

0 

(10c) 

3f 2 

30 ( r -°' 2 > 

= 0 

(10d) 


In a similar manner, conditions (7) when used in conjunction 
with (4), give rise to the following relations for 0 = tt , 
i . e . , 


3G 2 ( 1 - v ) 3f 3 . T o (r ’ z) 

3z r 30 2y 


(11a) 


36 2 ( 1 - v ) 9f l _ n 

3 r ~ r 39 


(lib) 


n ? 1 Hi 2 ( 1 -v ) 3f 2 ? 3 f 3 92f 0 9 2f 0 

(l-2v) 9r " r 3 Q 2v 9z + 9r ^ + 3 z z 


+ F ^ + z (f^- + t^)f 3 = 0 


(11c) 


where G can be found from 
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1 3 f Q 3 f -i 7 3 f 

6 - - ( 1 - 2V ) f 2 - - gg- + gg f g-Q-J 6 = IT 


( 12 ) 


Further simplification may be achieved by setting 


G = 0 , 0 = TT 


(13) 


On account of equation (13), it follows that 


3G _ 3_G 
3 z 3r 


0 , 0 


TT 


(14) 


and as an immediate consequence, equations (11a) and (11b) 
yield the crack-surface conditions 

3 f o rt (r,z) 

mr “ Ut v ) ■ 0 = * < 15 > 

a f , 

= 0, 0 = TT (16) 


Equations (15) and (16) together with (10b) and (10c) provide 
the necessary relations for the evaluation of f-j and f^. The 
remaining functions fg and f 2 are obtained from equations 
(13), (15), (16) and the fact that the loading is applied at 

the point r = a, z = a. This yields 

3f n 

(l-2v)f 2 - = 0 , e = TT (17) 

From equations (10a) and (lOd), it is not difficult to verify 

that 
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(18) 


37 [(1 " 2v > f 2 


3f ( 

3y 


-] = o, e = o 


The mixed conditions (17) and (18) provide a relation between 


the potentials f ^ 


3 f 


and 


0 


3y 


which is given in equation (30). 


The next step in the analysis is to derive a second relation 
between the same potentials. This may be accomplished by 
using equations (11c) and (10). Utilizing the identities 


3 r 


Jr 7 


3_ 
3x » 


3 z ' 


1 

r 30 


TP 


3_ 

3y 
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(19) 


equation (lie) can be transformed into 


2(l-v) 


3f. 

W 


- (l-2v) 


3 2 f . 


3f - 

air 


+ 2v 


3 f - 
3 ~z" 


3y 


= 0 , 0 = TT 


32f , 

3p“ 


3y3e 


At this point, it is expedient to add the term (y 


3 2 f 3 
3y3z 


) 


( 20 ) 


(which vanishes in the region 0 = tt ) to both sides of equa- 
tion (20). The resulting relation takes the form 


★ 

The second derivative of f^ is assumed to be regular for 

6 = m. This condition, incidently, is easily confirmed once 
the potential f^ is known. Refer to equation (26). 


- 10 - 
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ly ( 1 -v)f ^ - 



3f •, 3f o 

L + v — _ _ 

30 y 3 Z 



y 3fi 3fo 

- ( 1 - 2v) / (^ + 3 ^)dy] = 0, 9 = 7T (21) 

oo 


The above limits of integration have been introduced to ensure 
boundedness at the infinitely distant points. 

With reference to the region 9=0 and making use of some 
elementary transformations analogous to those in equations 
(19), equations (10) confirm that the relation (21) is also 
valid in the region 0=0. This means that equation (21) 
holds across the entire plane y = 0 of the half space y :> 0 . 

The expression inside the bracket in equation (21) is recog- 
nized as a harmonic function. It follows from Green's formu- 
la [6] that 

3 f n 3f, 3f o 3f, 

2 ( 1 -v ) f 2 - gy— = ge - + z 3^ y g"z 

y 3f, 3f o 

+ ( 1 - 2v ) / + ^)dy, y>0 (22) 


exists throughout the entire semi -infini te solid y^O. Equa- 
tion (22) provides the second equation for finding f 2 and f Q . 

METHOD OF SOLUTION 

In this section the potentials appearing in the basic 
displacement representation are determined by employing meth- 
ods of integral transforms and singular (primitive) solutions 
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of Laplace's equation [5]. 


The function f 3 (r,0,z) is governed by the mixed condi- 
tions (10c) and (15). This calls for the application of a 
Fourier cosine transform with respect to the variable z and 
a Kontorovi ch-Lebedev transform in r [4]. Seeking f 3 (r,e,z) 
by the expression 

f 3 = ( |) 7 / / *(s.t) K u (rs)cos(sz)dsdt (23) 


where is the Macdonald function [7] and 4>(s,t) is an ar- 
bitrary parameter such that equation (10c) is satisfied. By 
introducing the appropriate inversions and accounting for the 
concentrated load at (r = a, z = 0), 4> ( s , t ) is determined from 
equation (15) as 


<f>(s ,t) 


Rts i nh (irt ) K. t (rs ) 
( 2ir ) ^ ^tt 2 2y ( 1 - v ) 


(24) 


The result of putting equation (24) into (23) and employing 
the integral re presentation [7] is 


K. t ( rs ) K. t (as ) = / K Q ( s£ )cos ( tn ) dn 

(25) 

£ = (r 2 + 2arcos ho+a 2 )^^ 


* 


See Appendix for brief outline of this transform. 
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The integrals in equations (23) may now be evaluated to ren- 
der 


f 3 ( r , 0 , z ) = - 4 U7T ^( 1 . v ) p tan 1 (^ a ( r ~ x ) ) (26) 

where p represents the distance of any point in y^O from the 
point of application of R, i.e., 

p - :(x+a) 2 + y* + z ! ) 1/2 (27) 

The evaluation of f^(r,0,z) is dictated by the "homoge- 
neous" conditions (10b) and (16) governed by the Laplace 
equation in three dimensions. However, the choice f ^ = 0 is 
inappropriate as it does not lead to the correct solution 
(see [4], pages 381-383). This, in fact, is the crucial 
step in the present analysis and contradicts the general no- 
tion in the open literature that one of the harmonic func- 
tions contained in the Papkovi ch-Neuber solution can be 
dropped. In order to arrive at the proper expression for f ^ , 
a singular solution of Laplace equation which satisfies con- 
ditions (10b) and (16) and gives rise to well-behaved dis- 
placements at infinity is proposed [5] 

sin 4 

f,(r,6,z) = c -j Im[g(c)] (28) 

/r 

where Im designates the imaginary part of an analytic func- 
tion of the variable z defined by 
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?=r+a+iz, i = 


(29) 


and c-j is a constant introduced for convenience of writing 
the ensuing development. Note that as r->-0 , the expression 
(28) is undefined. On physical ground, however, the displace- 
ments must be finite in that region. This condition will be 
improved later on. 

To find the other potentials fg and f £, equations (17) 
and (18) suggest the following relation involving the special 
solution introduced in (28), i.e., 

0 

9fn A ^OS o' 

( 1 -2v)f 2 - j— = 2 Im[g ( c ) ] (30) 

y /r 


In equations (28) and (30) the function g(?) as well as the 
constants c-j and A are yet to be found. Solving equations 
(22) and (30) simultaneously and utilizing equations (26) and 
( 28 ) yield 


f 2 = { 4yfM-v) ' + 1 [ (3- 4 v)c-| - A]} 2 Im[g(?)] 

y sin 

+ (l-2v){i(»(x,y,z) - c, / Im[g'(e)]dy} (31a) 

oo /r 

and 


W = + 1 I(l-2v)(3-4v)c, - 2(1 -v)A]> 


cos 2 


Im[g(c)] + (1 -2v) (i|j ( x ,y , z ) - / 


y sin 


/F 


X 


Im[g ' (?) ]dy } 
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(31b) 



in which g'(?) = ^ g(?) and the following notations have 
been adopted: 


, . R k 0 l" vz -1 

* (x,y ’ z) = '4 'pttMI-v) tan 

- (2a) 1/2 Im[c o (2x-c o ) 
5 0 = ^ ] 9=0 = x + a + iz 


/ /2a ( r-x ) \ 

P ] 

i / o /r+x + /2x-r 

-V2 sn( 0 


) ]> (32 ) 


Integrating equation (31b) with respect to y between the 
limits °° and y and applying the results 


/ Im[g'(?)]dy = — Im j 

00 /r \/2 00 /t - c 

0 

1 '('(x.y.zjdy = n * Wl ; - / iU v) tf tan ’ 1 /MlEEI 

i/p 1 1/2 

- /Za Im[(c J-'^tan" 1 !— f-) ]} (33) 

^ ^0 

Q 

y y sin 7 . y ,,,,, 

/ dy / Im[g'(e)]dy = Im[ / 

00 00 vr 2/2 co /t-r 

^0 

- 2 y J g( t )^(t-? 0 + 2x)" 1/2 dt] 


it is found that 
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'o - - (^) ] 

/2y7T b s 0 


1/2 


+ [d-2v)c, - a] im ) 

/2 1 oo t-c 0 

+ + v ,Ml-v)p t 9 "' 1 (~ 2a - p r ' x) ) 


C 1 y ^ H 1/9 

+ — Im / 9 ( t ) jr( t-r +2x) 1/2 dt] 
/2 00 


(34) 


The remaining boundary condition to be satisfied is equation 
(10a). Imposing this condition on equation (34) results in 

[ ( 1 - 2v)c, - A] / SlUdt = _ ( . l“2v)R/a -1/2 ( 35 ) 

1 J r syrr* 2y ( 1 -v ) tt ‘’o ^ ; 

^ 0 ^0 

which is a standard integral equation of the Abel type [8] and 
thus the function g(?) is readily determined as 


g(?> - l 


(36) 


provided that 


(1 -2v)c-j - A = - 


(l-2v)R/a 

2y ( 1 - v ) tt z 


(37) 


The expressions for the potentials f Q and f 2 can be simplified 
further by defining a constant c 2 via the relation 


_ _ vR , r 1 \ 

C 2 “ 2yir 2 ( 1 - v ) ( 1 " v ) c i 


( 38 ) 


- 16 - 


Upon substituting equations (36-38) into (31a) and (34) the 
following expressions are derived 


Rz 


1 


fg( r ’6 »z) = (l-2v) [ 4y7r /ji _ v ) p tan ' x 


* (- 


2(ar) 1/2 sinf 


> + y4,] + ^ Im[ TviT 


Q/0 /r+x+/2x- F~ 

+ (2x-r )" 3/2 £n ( °)] 


(39) 


and 


f 2 (r,e,z) = c 2 


Q 

COS a“ ^ 

= Co Im(^r) + (l-2v){\f)(x,y,z) 


/r c 

Cl /r+x 


1 T_r fr+x _ \ -3/2 

- Im[ cT2x- ?0 ) - (2x ' 5 o) 

/r+x+/x- r 


x £n(- 


-)]> 


(40) 


where ip is defined in equation (32). The final step in the 
analysis is to obtain another expression relating c-| and c 2 - 
This may be done by imposing the regularity condition on the 
displacements as r->0. Inserting equations (26), (28), (39) 
and (40) in (2) and (3) and carrying out the expansion for 
small r, making use of (38) and retaining the lowest-order 
terms, it is found that the displacements near the crack edge 
become 
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(41a) 


fi 

c-.+c ? zsin p- 

u = . _J £ ■ / [3-4v + ( 7-8v ) cos 0 ] + 0(r°) 

2/r 

n 

c,+c ? zcos y 

U e = a *+ 7 * [3-4v - (5-8v)cose] + 0(r°) (41b) 

2/r 

u z = 0 ( r° ) (41c) 


The finite displacements requirement at the crack edge gives 


c 1 + c 2 = 0 (42) 

It follows from equation (38) that 
C 1 = " c 2 = " 2yTr r (1-vf(2-v) 

and the solution is complete. The physical quantities of in- 
terest may be readily computed from equations (2) - (5) when 
the appropriate expressions for the harmonic functions in equa- 
tions (26), (28), (39), (40) and (43) are used. 

STRESSES NEAR CRACK EDGE 


The shear stresses across the surface 0=0 are computed 
from equations (4b) and (4c) as 


T r0 (r ,° ,z ) 


4vR/a z (x+a ) 

^ (2 " v) (x ) 1 /2 [ (x+a ) 2 +z 2 ] 2 


T 0 z (r,o,z) 


R/a ( 2-3v) (x+a ) 2 + ( 2+v ) z 2 
irZ(2 " v) (x ) 1 /2 [ (x+a ) 2 + z 2 ] 2 


(44a) 

(44b) 
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Equations (44) may be expressed in the standard form as 


x r0 (r,o,z) 

T Qz (r,o,z) 


k 0 

— + 0 (x° ) 
/27 



0(x°) 


(45a) 

(45b) 


The stress-intensity factors for the edge-sliding and tearing 
modes of crack extension are expressed in terms of the non- 
dimensional parameter z* = through the relations 


k 


2 


4/2R , v x z* 

3/2 [ 2-v } (1+z**)* 
ua ' ' 


(46a) 


k 3 


/2R 


2 a 3/2 (l+z* 2 ) 


[1 


2v 


2- v 1+z 


l-z* 2 n 

irrl 


(46b) 


The variations of these equations with z* are shown in Figures 
2 and 3 for various values of v. 


CLOSURE 

The problem of determining stresses and displacements in 
a three-dimensional elastic solid with a semi -i nf i n i te plane 
crack is reduced to the finding of four distinct harmonic 
functions. The surfaces of the crack are deformed by the ap- 
plication of concentrated shears (R) parallel to the rectilin- 
ear boundary as shown in Figure 1. The analysis reveals that 
none of the Papkovi ch- Neuber potentials appearing in the gen- 
eral solution of the equations of equilibrium can be neglected. 
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Indeed, all of the faouji functions are needed to satisfy the 
appropriate boundary and regularity conditions in the problem. 

Methods of integral transforms involving the Macdonald 
function (the so-called Kontorovi ch-Lebedev transforms) as 
well as singular solutions of Laplace's equation in three di~ 
mensions are employed. The results presented here coupled 
with the solutions of Uflyand [4] complete the construction 
of the Green's functions for the semi-i nf i ni te crack problem 
and pave the way for immediate application in a number of 
problems of interest in fracture mechanics. 
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APPENDIX 


The Kontorovi ch-Lebedev transform of a function <f>(x) is 
defined by the formula 

A co , . 

4>(t) = / ^ K. t (x)dx (47) 

where <J) ( x ) is assumed to be of class L 2 in the interval [0,“] 

and K^(x) is the Macdonald function which is a solution of 

(’-£>» ■« <«> 


Assuming that 
•d> (x ) 


o(x) 


is continuously differentiable and x d> ( x ) 


and x — -] are absolutely integrable in [0,<»] then the 

inverse of (47) is given by 


- /s 

*(x) = ~~z / <|) (t ) K. . (x)tsinh (irt )dt 
IT O 1 Z 

Further details may be found in [4]. 


(49) 
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Figure 3. Variations of along crack border. 


